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We study dissipative transport of spontaneously emitting atoms in a ID standing-wave laser field 
in the regimes where the underlying deterministic Hamiltonian dynamics is regular and chaotic. 
A Monte Carlo stochastic wavefunction method is applied to simulate semiclassically the atomic 
dynamics with coupled internal and translational degrees of freedom. It is shown in numerical 
experiments and confirmed analytically that chaotic atomic transport can take the form either of 
ballistic motion or a random walking with specific statistical properties. The character of spatial and 
momentum diffusion in the ballistic atomic transport is shown to change abruptly in the atom-laser 
detuning regime where the Hamiltonian dynamics is irregular in the sense of dynamical chaos. We 
find a clear correlation between the behavior of the momentum diffusion coefficient and Hamiltonian 
chaos probability which is a manifestation of chaoticity of the fundamental atom-light interaction 
in the diffusive-like dissipative atomic transport. We propose to measure a linear extent of atomic 
clouds in a ID optical lattice and predict that, beginning with those values of the mean cloud's 
momentum for which the probability of Hamiltonian chaos is close to 1, the linear extent of the 
corresponding clouds should increase sharply. A sensitive dependence of statistical characteristics 
of dissipative transport on the values of the detuning allows to manipulate the atomic transport by 
changing the laser frequency. 
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I. INTRODUCTION 



An atom placed in a laser standing wave is acted upon 
by two radiation forces, deterministic dipole and stochas- 
tic dissipative ones [l|-y| • The mechanical action of light 
upon neutral atoms is at the heart of laser cooling, trap- 
ping, and Bosc-Einstcin condensation. Numerous appli- 
cations of the mechanical action of light include isotope 
separation, atomic litography and epitaxy, atomic-beam 
deflection and splitting, manipulating translational and 
internal atomic states, measurement of atomic positions, 
and many others. Atoms and ions in an optical lattice, 
formed by a laser standing wave, are perspective objects 
for implementation of quantum information processing 
and quantum computing. Advances in cooling and trap- 
ping of atoms, tailoring optical potentials of a desired 
form and dimension (including one-dimensional optical 
lattices) , controlling the level of dissipation and noise are 
now enabling the direct experiments with single atoms to 
study fundamental principles of quantum physics, quan- 
tum chaos, decoherence, and quantum-classical corre- 
spondence (for recent reviews on cold atoms in optical 
lattices see Ref. @, [![). 

Experimental study of quantum chaos has been car- 
ried out with ultracold atoms interacting in ^-kicked op- 
tical lattices - To suppress spontaneous emission 
(SE) and provide a coherent quantum dynamics atoms 
in those experiments were detuned far from the optical 
resonance. Adiabatic elimination of the excited state am- 
plitude leads to an effective Hamiltonian for the center- 
of-mass (CM) motion 14|, whose 3/2 degree-of- freedom 
classical analogue has a mixed phase space with regular 
islands embedded in a chaotic sea. De Brogile waves of 5- 



kicked ultracold atoms have been shown to demonstrate 
under appropriate conditions the effect of dynamical lo- 
calization in momentum distributions which means the 
quantum suppression of chaotic diffusion |5-[T3[. Deco- 
herence due to SE or noise tends to suppress this quan- 
tum effect and restore classical-like dynamics |15l - tl8| . 
Another important quantum chaotic phenomenon with 
cold atoms in far detuned optical lattices is a chaos as- 
sisted tunneling. In experiments (l9l - l2~l| ultracold atoms 
have been demonstrated to oscillate coherently between 
two regular regions in mixed phase space even though the 
classical transport between these regions is forbidden by 
a constant of motion (other than energy). 

The transport of cold atoms in optical lattices has been 
observed to take the form of ballistic motion, oscillations 
in wells of the optical potential, Brownian motion [22j], 
anomalous diffusion and Levy flights (23rl27| . The Levy 
flights have been found in the context of subrecoil laser 
cooling [23|, [24[ in the distributions of escape times for 
ultracold atoms trapped in the potential wells with mo- 
mentum states close to the dark state. In those experi- 
ments the variance and the mean time for atoms to leave 
the trap have been shown to be infinite. 

A new arena of quantum nonlinear dynamics with 
atoms in optical lattices is opened if we work near the op- 
tical resonance and take the dynamics of internal atomic 
states into account. A single atom in a standing-wave 
laser field may be treated as a nonlinear dynamical sys- 
tem with coupled internal (electronical) and external 
(mechanical) degrees of freedom (28rl30j . In the semi- 
classical and Hamiltonian limits (when one treats atoms 
as point-like particles and neglects SE and other losses of 
energy), a number of nonlinear dynamical effects have 
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been analytically and numerically demonstrated with 
this system: chaotic Rabi oscillations [28l - l30| . Hamil- 
tonian chaotic atomic transport and dynamical fractals 
j3ll - l33j , Levy flights and anomalous diffusion [30|, HH HH ■ 
These effects are caused by local instability of the CM 
motion in a laser field. A set of atomic trajectories un- 
der certain conditions becomes exponentially sensitive 
to small variations in initial quantum internal and clas- 
sical external states or/and in the control parameters, 
mainly, the atom-laser detuning. Hamiltonian evolution 
is a smooth process that is well described in a semiclassi- 
cal approximation by the coupled Hamilton-Schrodingcr 
equations. A detailed theory of Hamiltonian chaotic 
transport of atoms in a laser standing wave has been de- 
veloped in our recent paper [33| . 

The aim of the present paper is to study dissipative 
chaotic transport of atoms in a one-dimensional optical 
lattice in the presence of SE events which interrupt co- 
herent Hamiltonian evolution at random time instants. 
Generally speaking, deterministic (dynamical) chaos is 
practically indistinguishable in some manifestations from 
a random (stochastic) process. The problem becomes 
much more complicated when noise acts on a dynamical 
system which is chaotic in the absence of noise. Such 
systems are of a great practical interest. Comparatively 
weak noise may be treated as a small perturbation to 
deterministic equations of motion, and one can study 
in which way the noise modifies deterministic evolution 
on different time scales. However, SE is a specific shot 
quantum noise that cannot be treated as a weak one be- 
cause the internal state may change significantly after the 
emission of a spontaneous photon. Special methods are 
needed to describe correctly the dynamics of a sponta- 
neously emitting single atom in an optical lattice. The 
purpose of this paper is twofold. Our first goal is to give 
a description of possible regimes of dissipative atomic 
transport in the presence of SE and to quantify their sta- 
tistical properties. Our secondary intent is a search for 
manifestations of the fundamental dynamical instability 
and Hamiltonian atomic chaos in the diffusive-like CM 
motion of spontaneously emitting atoms in a laser stand- 
ing wave which can be observed in real experiments. 

The paper is organized as follows. In Sec. II we for- 
mulate a Monte Carlo stochastic wavefunction approach 
to solving semiclassical Hamilton-Schrodinger equations 
of motion for a two-level atom in a one-dimensional 
monochromatic standing light wave. This approach al- 
lows to get the most probabilistic outcome that can be 
compared directly with corresponding experimental out- 
put with single atoms. In Sec. Ill we review briefly 
our previous results on Hamiltonian chaotic CM motion 
which are necessary to quantify and interpret dissipative 
dynamics. Sec. IV is devoted to description of possi- 
ble regimes of dissipative CM motion of spontaneously 
emitting atoms in a standing wave. Monte Carlo sim- 
ulations of the well-known effects of acceleration, decel- 
eration, and velocity grouping, and of a novel effect of 
dissipative chaotic walking of atoms are presented in this 



section. Anomalous statistical properties of dissipative 
chaotic walking are quantified and discussed in Sec. V. 
Whereas Sees. IV and V are devoted to solving the first 
task of this paper, in Sec. VI we consider the problem of 
manifestations of dynamical instability and Hamiltonian 
chaos in dissipative atomic transport. We demonstrate 
analytically and numerically that character of diffusion 
of spontaneously emitting atoms changes qualitatively in 
the detuning regime where the underlying Hamiltonian 
dynamics is chaotic. To observe this effect in a real ex- 
periment with cold atoms in a one-dimensional optical 
lattice we propose to measure the linear extent of atomic 
clouds with different values of their mean momentum and 
predict that the extent should increase significantly with 
those values of the mean momentum for which the un- 
derlying Hamiltonian evolution is chaotic. 



II. MONTE CARLO WAVEFUNCTION 
MODELING OF THE ATOMIC DYNAMICS 



In the frame rotating with the laser frequency w/, the 
standard Hamiltonian of a two-level atom in a strong 
standing-wave ID laser field has the form 
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where <t± i2 are the Pauli operators for the internal atomic 
degrees of freedom, X and P are the atomic position and 
momentum operators, uj a , uif, and Q arc the atomic tran- 
sition, laser, and Rabi frequencies, respectively, and T is 
the spontaneous decay rate. Internal atomic states are 
described by the wavefunction \^{t)) = a(t)\2) + 6(t)|l), 
with a and b being the complex-valued probability am- 
plitudes to find an atom in the excited \2) and ground |1) 
states. Note that the norm of the wavefunction, |a| 2 + |6| 2 , 
is not conserved due to non-Hermitean term in the Hamil- 
tonian. 

We use the standard Monte Carlo wavefunction tech- 
nique [36[ to simulate the atomic dynamics with the cou- 
pled internal and external degrees of freedom in an opti- 
cal lattice. The evolution of an atomic state |\P(t)) con- 
sists of two parts: (i) jumps to the ground state (a = 0, 
\b\ 2 = 1) each of which is accompanied by the emission 
of an observable photon at random time moments with 
the mean time 2/T (actually, the probability of SE de- 
pends on the atomic population inversion) and (ii) co- 
herent evolution with continuously decaying norm of the 
atomic state vector without the emission of an observ- 
able photon. The decaying norm of the state vector is 
equal to the probability of spontaneous emission of the 
next photon. It is convenient to introduce the new real- 
valued variables (normalized all the time) instead of the 
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amplitudes a and b (rcnormalizcd after SE events only) 
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(2) 

which have the meaning of synphase and quadrature 
components of the atomic electric dipolc moment and 
the population inversion, respectively. We stress that the 
length of the Bloch vector, u 2 + v 2 + z 2 = 1, is conserved. 

Since we study manifestation of quantum nonlinear ef- 
fects in ballistic transport of atoms, when the average 
atomic momentum is very large as compared with the 
photon momentum hkf, the translational motion is de- 
scribed classically by Hamilton equations. The whole 
atomic dynamics is governed by the following Hamilton- 
Schrodinger equations [ID, |43| 

oo 
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where x = kf(X) and p = (P)/Hkf are normalized 
atomic CM position and momentum. The dot denotes 
differentiation with respect to the normalized time r = 
fit. Throughout the paper we fix the values of the nor- 
malized decay rate 7 = T/fi and the recoil frequency 
ui r = hk 2 f /m a n to be 7 = 3.3 ■ 1CT 3 and ui r = 1CT 5 . This 
values are similar to those used in experiments with Na 
@, 0|, Cs 0, S3 and Rb [H] cold atoms in a standing- 
wave laser field with the maximal Rabi frequency of the 
order of 1-^5 GHz. So, the normalized detuning between 
the field and atomic frequencies, A = (ojf — oj a )/Q, is a 
single variable parameter. Also we fix the initial condi- 
tions as follows: xq = vq = uq = 0, zq = —1, and vary 
the initial momentum po only. In Eqs. §3§ Tj are random 
time moments of SE events and pj are random recoil mo- 
menta with the values between ±1 (ID case). In terms of 
the normalized time t the rate of occurrence of SE events 
is 7(2; + l)/2. At t = Tj, the atomic variables change as 
follows: 



0, v — > 0, z — > —1, p — > p+pj, 



T <Pj < 1. 

(4) 



III. 



A BRIEF REVIEW OF HAMILTONIAN 
ATOMIC DYNAMICS 



In this section we review briefly our main results on 
Hamiltonian atomic dynamics (sec Rcfs. 28 31, 33, 35]) 



which will be used in the next sections. In the absence of 
any losses (7 = 0) the total atomic energy is conserved, 
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The corresponding lossless equations of motion with two 
independent integrals of motion, the energy H and the 
length of the Bloch vector, have been shown [2^, [29| to 
be chaotic in the sense of an exponential sensitivity to 
small variations in initial conditions and/or the control 
parameters. The CM motion is governed by the simple 
equation for a nonlinear physical pendulum with the fre- 
quency modulation [39| 



x + cj. r u(r) sinx = 0, 
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where the synchronized component of the atomic dipolc 
u is a function of all the other atomic variables includ- 
ing the translational ones. Besides the regular CM mo- 
tion, namely, oscillations in a well of the optical potential 
and a ballistic motion over its hills, we have found and 
quantified chaotic CM motion [2^, [2i| 39] . On the exact 
atom-laser resonance with A = 0, u is a constant, and 
the CM performs either regular oscillations, if H < \u\, 
or moves ballistically, if H > \u\. 

At A ^ 0, the depth of the potential wells changes in 
course of time, and atoms may wander in a rigid optical 
lattice (without any modulations of its parameters) in a 
chaotic way with alternating trappings in the wells and 
flights of different lengths and directions over the hills. 
At small detunings, |A| <C 1, the second term of the 
potential energy U in Eq. ([5]) is small, and U can be 
approximated by a function of only one internal variable 
u. In this case we have approximate solutions for v and 
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where xo is an integration constant which is a function 
of initial values of z and u. Using these solutions one can 
prove that at |A| <C 1, u performs shallow oscillations 
when the atom moves between the nodes (recall that u = 
const at A = 0). These oscillations are synchronized with 
the oscillations of z, and when an atom approaches any 
node with cos a; = 0, where the strength of the laser field 
changes the sign, they slow down (see Eq. (J7J). The 
swing of oscillations of u gradually increases, and exactly 
at the node u changes abruptly its value (see Fig. [T]). 
Thus, u is practically a constant between the nodes and 
it performs a sudden jump at every node. 

In the Raman-Nath approximation, where x = tu r pT 
and p = const, we have managed to derive the deter- 
ministic mapping allowing to compute the value u m just 
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after crossing the mth node 
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z at the antin- 
k = 0,1,2,.... 



where vq and z$ are the values of v and 
odes of the standing wave at x = irk, 
They are the same at all the antinodes because in the 
Raman-Nath approximation v and z are periodic func- 
tions of x (see solution ([?]))• Formula © describes the 
series of jumps of two alternating magnitudes (for odd 
and even m). Strictly speaking, (|8]) is valid with fast bal- 
listic atoms and not on a very long time scale. Deviation 
of the analytic calculations with Eq. from the exact 
numerical results is demonstrated in Fig. QJi where we 
plot the function u(t) for a fast atom with po = 1900. It 
is obvious that the signal is rather regular but the mag- 
nitude of the jumps changes slowly because the Bloch 
components v and z are not strictly periodic functions of 
time. 

Figure [TJd plots u(r) in the regime of Hamiltonian 
chaotic walking. To quantify chaotic jumps of u we pro- 
posed in Ref. [33| the stochastic map 
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which was derived from the deterministic map (|5]) by in- 
troducing random phases <p m (0 < (/) < 27r) instead of 
arguments of trigonometric functions which may differ 
significantly from node to node due to strong variations in 
the atomic momentum beyond the Raman-Nath approx- 
imation. Note that the value of the momentum at the 
instant when the atom crosses a node, p no dc = \/2H/uj r , 
is approximately the same for all nodes. 

The map describes a random Markov process in 
the u space with u m varying in the range — 1 < u m < 1. 
This quantity may be smaller or larger than the atomic 
energy H (which is a constant in the Hamiltonian limit). 
Since the values of u define the atomic potential energy, 
its random walking governs a random walking of atoms in 
the lattice. The possible regimes of the Hamiltonian CM 
motion can be summarized as follows [33| : At \u\ > H, an 
atom oscillates in one of the potential wells, at |u| < H . it 
moves ballistically. It can walk chaotically if < H < 1. 
In the process of Hamiltonian chaotic walking the atom 
wanders in an optical lattice with alternating trappings 
in wells of the optical potential and flights over its hills 
changing the direction of motion many times. "A flight" 
is an event when the atom passes, at least, three nodes. 
CM oscillations in a well of the optical potential is called 
"a trapping" . The number of node crossings I during a 
single flight or a single trapping event is a discrete mea- 
sure of the length and durations of those events. We 
have derived in Ref. [H| the following formulas for the 



probability density functions (PDFs) for the flight and 
trapping events in the diffusive approximation: 
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Here Q is a constant, D u = A 2 7r/4w r .p no do is a diffusion 
coefficient in the u space. For comparatively small values 
of I (i. e., with short flights and trappings), we get from 
Eq. (fl~U)) the power decay 
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whereas for large I the decay is exponential. Numerical 
simulation of the Hamiltonian equations of motion agrees 
well with the analytical results (pl7|) in different ranges of 
the detunings. A typical PDF for the flight and trapping 
events decays initially algebraically and has an exponen- 
tial tail. The length of the initial power-law segment is 
inversely proportional to the value of the detuning A and 
can be rather large. 

In which way SE changes the statistical properties of 
the Hamiltonian motion? Can we find fingerprints of 
Hamiltonian instability and chaos in the motion of spon- 
taneously emitting atoms or SE totally suppresses any 
manifestations of coherent (but chaotic!) Hamiltonian 
dynamics? These questions will be addressed in the next 
sections. 



IV. DISSIPATIVE ATOMIC TRANSPORT IN A 
LASER STANDING WAVE 

The emission of a photon into the continuum of modes 
of the electromagnetic field is accompanied by an atomic 
recoil. The dissipative (friction) force F = (p) (which 
does not exist in the Hamiltonian system) depends on 
the atomic momentum p and the sign of the detuning in 
a complicated way 0, |4l|. The effects of acceleration, 
deceleration, and velocity grouping (at A < 0) are well- 
known in the literature [2, 3]. A novel effect we report in 
this section is dissipative chaotic walking. It appears un- 
der appropriate conditions that arc different from those 
specified for Hamiltonian chaotic walking in the preced- 
ing section. 

To illustrate the possible regimes of dissipative atomic 
transport in a standing wave we integrate by the Monte 
Carlo method dissipative equations of motion ([3]) with 
2000 atoms whose positions and momenta arc distributed 
in a quasi-Gaussian manner (Fig. H£i). In Fig. ^Bp we 
demonstrate the velocity grouping effect at A = —0.2 and 
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FIG. 1. Time evolution (r is in units of of the synphase component of the electric dipole moment u and the atomic 

energy H (in units of ftfl). (a) Regular Hamiltonian dynamics (po = 1900, 7 = 0), (b) chaotic Hamiltonian dynamics (po = 700, 
7 = 0), (c) regular dissipative dynamics (po = 1900, 7 = 0.0033), (d) chaotic dissipative dynamics (po = 700, 7 = 0.0033). In 
all the panels, A = —0.0005. The initial part of (a) agrees with approximate solution ([8]) with vq = 0, zq = —1. 



t = 10 5 . A large number of atoms is grouped around two 
values of the capture momentum p cap ~ ±1300 because 
of acceleration of slow atoms and deceleration of the fast 
ones in the initial ensemble. The slower the atoms are the 
longer is the process of the velocity grouping. Note that 
atoms with \p\ < 100, trapped initially in a well of the 
optical potential, could not quit the well up to r = 10 5 . 
Contrary to that, at positive values of the detuning fast 



atoms are accelerated and slow ones are decelerated. As 
a result, we observe a pronounced peak around x ~ p ~ 
shown in Fig. at A = 0.1 and r = 10 5 . 

Dependence of the friction force F on the current 
atomic momentum p is shown in Fig. [3] at A = —0.2. It 
has been computed with our main equations © when av- 
eraging over seven thousands atoms with different initial 
momentum. The function F(p) resembles the behavior 
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FIG. 2. Atomic momentum and position distribution illustrating the effects of atomic acceleration, deceleration, and the 
velocity grouping: (a) r = 0, (b) r = 10 s , A = -0.2, (c) r = 10 5 , A = 0.1, (d) r = 10 5 , A = -0.05. Momentum p is given in 
units of Hkf, the position in units of kj 1 . 



of the friction force computed with another methods (see 
and Fig. la in [4lj). The friction force decreases up to 
zero value and then begin to oscillate with increasing p. 
It has a number of zeroes (the detailed view is shown in 
Fig-Eb) like the corresponding functions in Refs. [3, l4lj. 
Zero values of F correspond to quasistationary values of 
the momentum which depend on A. Some of them are 
attractors and atoms with close values of the momentum 
tend to Pcap; another ones arc rcpcllors. The attractors 
and rcpcllors arc not deterministic because of a random 
nature of SE. Thus, an atoms walks randomly in the mo- 
mentum space between different values of the capture 
momentum p cap . When it reaches a certain value of the 
capture momentum the atom does not stop in the mo- 
mentum space and goes on to fluctuate because of both 
the Hamiltonian instability and SE effect. 

In the preceding section we described the Hamiltonian 
chaotic walking that may occur in the absence of any 
losses. Dissipation causes additional strong fluctuations 
of the momentum. If A > or if it is negative but 
comparatively large, nothing principally new happens to 
atoms as compared with the Hamiltonian limit. How- 
ever, at negative small values of A, a characteristic cap- 
ture momentum becomes smaller than a typical range of 
momentum fluctuations due to atomic recoils. As a re- 



sult, atoms may change their direction of motion in an 
irregular way. Such a dissipative chaotic atomic walking 
is illustrated in Fig. EJl at A = -0.05 and r = 10 5 with 
the atoms distributed widely in the phase plane. Typical 
atomic trajectories are shown in Fig.Uin the momentum 
space. Figures H^, and b illustrate how the friction force 
near the resonance (A = —0.001 and A = —0.01) decel- 
erates atoms with large values of the initial momentum 
down to so small values of the capture momentum when 
the dissipative chaotic walking becomes possible. With 
increasing the absolute value of the negative A, the cap- 
ture momentum increases and the atom changes rarely 
the direction of motion (Fig. 0]; with A = —0.1). Panels 
d and e in Fig. [4] illustrate the velocity grouping effect at 
A = —0.15 with different values of the initial momentum. 



V. STATISTICAL PROPERTIES OF 
DISSIPATIVE CHAOTIC WALKING 

Statistics of Hamiltonian chaotic walking is quanti- 
fied by the flight and trapping PDFs (fT0|) with alge- 
braically decaying head segments and exponential tails 
whose lengths strongly depend on A. We will show in 
this section that PDFs for dissipative chaotic walking 
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FIG. 3. Dependence of the friction force F = (p) on the 
current atomic momentum p at the detuning A = —0.2. 



are even more sensitive to variations in A. Figures |4^, 
and b clearly demonstrate that at very small detuning 
A = —0.001 long flights dominate, whereas there appears 
a number of short flights with larger value of A = —0.01. 
In Fig. [5]we plot the PDFs P/i for the duration of atomic 
flights T with different values of the detuning A. At 
small dctunings (Fig. [5^,), the length of the power-law 
segments depends on A in a similar way as in the Hamil- 
tonian case (compare this figure with Figs. 5a, 6a, and 
7a in Rcf. [33]). However, the slope slightly differs from 
the Hamiltonian slope which is equal to —1.5. The dif- 
ference in the statistics of dissipativc and Hamiltonian 
walkings is more evident with larger values of the de- 
tuning (Fig. [SJd). The length of the power-law segments 
increases drastically with increasing A. This effect is ab- 
sent in Hamiltonian dynamics. The corresponding slope 
a decreases with changing the detuning from A = —0.09 
to A = —0.12 because of the corresponding increase in 
the length of atomic flights (see Figs. [4^, b, and c). In 
Fig. |5] we plot the dependencies of the mean duration 
of atomic flights (T) and the slope of the PDF power- 
law fragments a on the detuning A in the range of its 
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FIG. 4. Typical atomic trajectories in the momentum space: 
(a)-(c) dissipative chaotic walking with different statistics of 
atomic flights, (d)-(e) the effect of velocity grouping. Note 
that atoms with very different initial momentum acquire a 
close value of the capture momentum. 



medium values —0.12 < A < —0.06. Both the quantities 
correlate well with each other. It means that, changing 
the value of A, one can manipulate statistical properties 
of dissipative atomic transport in an optical lattice. This 
control is nonlinear in the sense that slightly changing A, 
say, from —0.08 to —0.12, we increase the mean duration 
of flights in a few orders of magnitude. This effect may 



8 



8# X 



CL 

o 

o 



-10 



-12 




-12 




4 5 6 7 8 

iogi T 

FIG. 5. The PDFs P fl for the duration of atomic flights T with 
(a) small detunings (crosses A = —0.01, stars A = —0.001, 
circles A = -0.0001, squares A = -0.00001) and (b) medium 
detunings (stars A = -0.09, a = -0.77; circles A = -0.1, 
a = -0.27; squares A = -0.12, a = -0.05). Straight lines 
show slopes a of the power-law fragments of the PDFs in 
log-log scale. 



be qualitatively explained as follows. When increasing 
the absolute value of the negative detuning, the capture 
momentum increases but fluctuations of the current mo- 
mentum p decrease providing long atomic flights [43j . 

To explain the statistical properties of the dissipa- 
tive chaotic walking let us consider the behavior of the 
quasienergy 



Hi 



UJ r 2 A A7 2 

—p -UCOSX-—Z-—{l-Z){T-T j ) = 

= H-^-{l-Z 2 )(T-Tj), Tj < T < Tj—i , 

(12) 

which is equal to the total atomic energy ((SJ) in the ab- 
sence of relaxation. Near the resonance, |A| <C 1, Hj 
is almost conserved between SE events, i. e., in the in- 
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FIG. 6. Dependencies of the logarithms of the mean duration 
of atomic flights (T) (solid line) and of the slope a of the PDF 
power-law fragments (squares) on the detuning A (in units of 
Q). 



terval tj < t < Tj-!. The real energy H (see Fig. Qji) 
decreases a little in between in a linear way. The rate of 
this decrease is defined by the coefficients of spontaneous 
emission 7, the detuning A, and the average probability 
to find the atom in the excited state. Both the quantities, 
H and H, changes abruptly just after SE (because of the 
corresponding changes in the atomic variables ((3])). Just 
after emitting a jth spontaneous photon at r = Tj , they 
have the same values. So, we will model the evolution of 
the energy as a map Hj = H{t^) taken at the moments 
rf just after SE 



Hj = Hj - Hj-! + Hj-! = Hj-! + u r p(Tj ) Pj + 

U r n A 

+ Y p 3 + ~2 + U ( T i ' C0SX ( T j) + 

+ J Z (Ti) + ^ L (l-Z 2 )(Tj-Tj- 1 ), 



(13) 



where the values of the atomic variables p(r- ),u{t- ), 
and z(t~) are taken at the moments tJ just before SE. 
They are, in turn, determined by the coherent evolution 
between SE. 

The stochastic map for the atomic energy (|13p provides 
an important information about the CM motion. It has 
been shown in Ref. [HJ that atoms move ballistically, 
if the atomic energy satisfies to the condition H > \u\, 
whereas at H < |u| they may change the direction of mo- 
tion. The dissipative chaotic walking takes place when 
the atomic energy H alternatively takes the values larger 
and smaller than a critical value H = \u\. In the Hamil- 
tonian limit, where the energy is conserved, the problem 
of the CM chaotic walking has been reduced to the task 
of random walking of the Bloch component u (see Sec. Ill 
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and Ref. [3J]). The energy is not conserved in the pres- 
ence of relaxation, but the values of u are always small 
(see Appendix and Fig. [TJ; and d). Thus, atoms oscillate 
in the wells of the optical potential if H < and move 
ballistically if H > 0. 

On a time scale exceeding the mean time between SE 
events 2/7, the evolution of energy can be treated as a 
diffusion process with a drift in the energetic space. The 
probability to have the energy H at time r is governed 
by the Fokker-Planck equation 



P(H,r) 



-2c H 



dP 
dH 



D 



H 



d 2 p 

dH 2 ' 



(14) 



where Dh is an energy diffusion coefficient and ch is an 
energy drift coefficient which can be estimated with the 
help of Eq. (p~3|) as follows: 



ch 



{Hi - g 3 -_!) 



(H) 



12 



A7 



(15) 



In deriving this formula we adopt the average value of 
the squared recoil momentum = 1/3 (the projec- 
tions of the recoil momenta on the standing-wave axis x, 
Pj, are assumed to be distributed in the range ±1 with 
the same probability), the average value of the popula- 
tion inversion just before a SE event to be z(t~) = 1/2, 
the average value of z over the whole time scale and 
its mean squared deviation from 1 to be (z) — and 
(1 — z 2 ) = 1/2, respectively (see the solution (f3"3")l in Ap- 
pendix). Moreover, neglecting the correlation, we put 
(ucosx) ~ (u)(cosx) = 0, which is valid if \ui r p\ > 7/2, 
i. e., when p > 100 with our choice of the parameters. 
Since the first term in (|15[) is small and may be neglected, 
the drift velocity of an atom in the energetic space is ap- 
proximately proportional to the detuning A, and, there- 
fore in average atoms accelerate and decelerate at A > 
and A < 0, respectively, as it should be for |A| <C 1. In 
the weak Raman-Nath approximation, (|32|) and (f34|) . the 
drift coefficient in the energetic space is simply related to 
the friction force F acting upon atoms 



F=(p) 



ui r p 



(16) 



The friction force plays the role of a drift coefficient in the 
momentum space. Strictly speaking, the weak Raman- 
Nath approximation is not valid near the turning points 
when the atomic velocity is comparatively small. How- 
ever, most of the flight time it is valid. 

The diffusion coefficient in the energetic space is given 
by the formula 



D 



H 



((Hj - Hj-1) 2 ) - (Hj - Hj-i)' 
2<r i -r J -_i> 



(17) 



which can be rewritten with the help of (|X3|) as follows: 
Dh ~ ^Vfcr) + (^(rr)h (18) 



Using weak Raman-Nath approximation, (|32|) and (|34[) . 
the first term can be replaced by juj r H/6. Using the 
estimation (|38j) for (u 2 {t~)) in the irregular CM motion 
regime (see Appendix), we get the following expression 
for the energy diffusion coefficient: 



D 



H 



^uj r H A 2 
6 + T 



(19) 



This expression is valid for moderately small momcntums 
(p 1000) when the strong Raman-Nath approximation 
cannot be applied. In the process of dissipative chaotic 
walking, the probability to get higher values of the mo- 
mentums is almost zero. 

Now we will try to derive analytically a distribution 
of the durations T of atomic flights in the process of 
dissipative chaotic walking. In fact, it is a problem of 
the first passage time for the atomic energy H to return 
to its zero value. Recall that at small dctunings we have 
H ~ in the very beginning of every flight. In course of 
time H can reach rather large values, and it returns to 
zero at the end of the flight. If the random jumps of the 
energy would be symmetric (ch = 0), the probability 
to find a flight with duration T would be proportional 
to T~ 1,5 , where the exponent —1.5 does not depend on 
the diffusion coefficient. This conjecture follows from the 
known theorem in probability theory. More general result 
(see chapter XIV in Ref. j44j) proves that in the case 
of an asymmetric random walking in the energetic space 
(eff 7^ 0) the PDF for the flight durations in configuration 
space is 



Pfl oc e -<? B T/D aT -i.s 



(20) 
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if the drift and diffusion coefficients in the Fokker-Planck 
equation for the random walking are assumed to be con- 
stants. This formula gives a distribution of the flight 
durations with a power-law fragment followed by an ex- 
ponential tail and agrees qualitatively with the exact nu- 
merical computations of Pfi shown in Fig. [S^i for a few 
values of the detuning A. The main disadvantage of this 
formula is that (f20|) does not depend on A as the ex- 
act PDFs in Fig. [SJi. At very small A = -10 -5 , the 
PDF, shown by squares in Fig. [5^, decays mostly alge- 
braically, whereas at larger values of the detuning the 
power-law fragments are much shorter. A discrepancy 
between the analytical and numerical PDFs arises be- 
cause we assumed in deriving (|20[) that Dh and ch do 
not depend on the energy H . In fact, it is not the case 
for small values of the momentum p, and a more accurate 
formula for Pfi(T) is required. 

The PDFs for Hamiltonian ((TUJ) and dissipative (|20|) 
transport are similar in the sense that both Pfi contain 
power-law fragments followed by exponential tails, but 
the origin of each statistics is different. In the Hamilto- 
nian limit the statistics is governed by the behavior of u, 
not the energy, as in the dissipative system. A turnover 
from a power law to an exponential decay in the Hamil- 
tonian case is explained by a boundedness of |u| < 1, 
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whereas in the dissipative system it is explained by a 
negative drift of the energy H. Each of the factors pre- 
vents the corresponding randomly walking quantity to go 
far away from its critical value (at which the atoms can 
change the directions of motion) decreasing the probabil- 
ity of long flights in at exponential way. 



VI. MANIFESTATION OF HAMILTONIAN 
CHAOS IN DISSIPATIVE ATOMIC TRANSPORT 

In the absence of SE the atomic dynamics can be regu- 
lar or chaotic depending on the initial conditions and/or 
the detuning. In experiments one measures statistical 
characteristics of spontaneously emitting atoms. Is there 
a correlation between those characteristics and the un- 
derlying Hamiltonian dynamics? Can we find any man- 
ifestations of Hamiltonian instability, chaos, and order, 
in the diffusive-like dissipative atomic transport? These 
questions will be addressed in the present section. 

The common quantitative criterion of deterministic 
chaos, the maximal Lyapunov exponent A, is a measure of 
a divergence of two trajectories in the phase space with 
close initial conditions [45} . To quantify probability of 
chaos in the mixed Hamiltonian dynamics, when A = 
with some values of po and A > with another values of 
Po, we introduce a probabilistic measure of Hamiltonian 
chaos 



A = (29(A) -1), 



(21) 



where 0(A) is a Heaviside function (6 = if A < 0, 
9 = 1/2 if A = 0, and 9 = 1 if A > 0). The probability 
of Hamiltonian chaos A is computed by averaging over a 
large number of atomic trajectories with different values 
of p . If all the trajectories in the set turn out to be 
stable, one gets A = 0, and if all they are exponentially 
unstable, then A = 1. One gets < A < 1, if some 
trajectories in the set are stable but the other ones are 
not. The magnitude of A is proportional to the fraction 
of trajectories with positive As. 

To examine manifestations of the underlying Hamilto- 
nian dynamics in dissipative transport in is convenient to 
consider atomic diffusion not in the energetic but in the 
momentum space. The momentum diffusion coefficient, 
which is a measure of momentum fluctuations, can be 
written with the help of (|17|) and (|52"j) as follows: 



D H 7 , ( u ( T j )>7 



ui 2 p 2 12 8w 2 p 2 



(22) 



Using the formula (fl9| , we get D p in the regime of chaotic 
oscillations of the Bloch component u 



,7 A 2 



12 8lu 2 p 2 ' 



(23) 



The momentum diffusion coefficient D p is computed with 
the main equations ([3]) in the following way. The range 



of all possible values of the atomic energy H §5§ is parti- 
tioned in a large number of bins. For a large number of 
initial conditions (in fact, we change only the initial mo- 
mentum po keeping the other conditions to be fixed), af- 
ter any jth SE event we compute the difference Hj—Hj—i 
and the squared difference {Hj — Hj-i) 2 . They are ran- 
dom values, but their statistics depend on the preceding 
energy value Hj-i. So we calculate the histograms of 
the average values of (Hj — ffj-i) and {{Hj — Hj-i) 2 ) 
as functions of energy Hj-\. After that we can compute 
the energy diffusion coefficient Dh (| 1 T[) which, being di- 
vided by u 2 p 2 , yields the momentum diffusion coefficient 
Dp which is better to present as a function of the current 
momentum p ~ ^/2Hj^i/uj r . 

The main result in this section is illustrated with Fig. [7] 
In the upper left and right panels the dependencies of the 
momentum diffusion coefficient D p on the current mo- 
mentum p are plotted in a log-log scale for A = —0.01 
and A = —0.0005, respectively. In both the cases, we 
put 7 = 0.0033. These plots should be compared with 
the corresponding lower panels where the probability of 
the Hamiltonian chaos A is plotted against p with 7 = 
(i. e. in the Hamiltonian limit of Eqs. <j3j) ) . It is evident 
that the character of the momentum diffusion changes 
abruptly at those values of the current momentum where 
a transition from chaos to order occurs in the underlying 
Hamiltonian dynamics. Such a turnover takes place in a 
range of small negative detunings and is a manifestation 
of the peculiarities of the underlying Hamiltonian evo- 
lution in the diffusive-like dissipative transport of atoms 
in a standing-wave laser field. We may conclude that in 
spite of random atomic recoils due to SE the chaotic (reg- 
ular) dynamics between the acts of SE clearly manifests 
itself in the behavior of the measurable characteristic of 
the atomic transport, the momentum diffusion coefficient 
Dp. The behavior of D p in the range of p, where the 
underlying Hamiltonian evolution is chaotic, is well de- 
scribed by the formula (|23|) with Dp 1 ~ p 2 + const (see 
both the upper panels in Fig. [7] where this dependence is 
shown by dashed lines). 

However, the formula (|23[) does not work in the regimes 
when the underlying Hamiltonian dynamics is mixed or 
regular because in deriving it we supposed fully chaotic 
behavior of u. We have managed to estimate analytically 
Dp in the Hamiltonian regular regime at extremely small 
values of the detuning | A | <^ 1 and for atoms whose mo- 
mentum is so large that we can neglect its fluctuations 
between SE events (the exact Raman-Nath approxima- 
tion with x = u r pr). Figure UK illustrates the ladder- like 
behavior of u which is descibed by the deterministic map- 
ping © on a comparatively short timescale. To get D^ N 
from Eq. (f2"2"j) we use the expression (f3"T| for u 2 {tJ ) de- 
rived in Appendix 



D* N ~^ 



1 



12 SuJrP'yir 



(24) 



Thus, we derived the formulas for the momentum dif- 
fusion coefficient D p in the regimes of Hamiltonian chaos 
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FIG. 7. Correlation between the average momentum diffusion coefficient in a log-log scale D p (in units of fi 2 fc?f2) and probability 
of Hamiltonian chaos A in their dependencies on current atomic momentum p (in units of hkf) at A = —0.01 and A = —0.0005. 
Dashed lines with the slopes p _2 +const and p _1 + const are theoretical curves ([23} and (|24[) valid in the regimes of Hamiltonian 
chaos and order, respectively. Solid line is a theoretical curve (|25J1 derived to fit the exact numerical results. An abrupt change 
in the decay laws for D p occurs for those values of p where the transition from order to chaos takes place in the underlying 
Hamiltonian dynamics. 



(123)) with A = 1 and Hamiltonian order ([24]) with A = 0. 
In a general case with < A < 1, we will assume a linear 
combination 



D p ~(l- A)D, 



RN 
P 



12 



A 2 



8u> r p \ 77T 



1 - A 



A 

LO r p 



(25) 



This function, shown by the solid line in the right upper 
panel in Fig. [71 fits rather well exact numerical results. 

We would like to end this section with the proposal of 
a simple experimental scheme to observe our main theo- 
retical and numerical result on an abrupt change in the 
character of atomic diffusion in a laser standing wave un- 
der conditions corresponding to two different regimes of 
the underlying Hamiltonian evolution, chaotic and reg- 
ular ones. Let us consider a small atomic cloud mov- 
ing in one direction with an average atomic momentum 
(p c )- Initial position and momentum distribution are as- 
sumed to be a Gaussian with the standard deviations 
a 2 = ((x — (x c )) 2 } and a 2 = ((p - (p c )) 2 ). The momen- 
tum diffusion coefficient is 



d(al) 



D *= 2dr 



(26) 



The temperature of gas and its rate of heating in Kelvins 



per second are 

rp _ 2{E k ) 



h 2 k}a 2 p 
m a k B 



dT 
~~dl 



2h 2 k 2 flD p 
m a k B 



(27) 



where Ek is a kinetic energy of atoms (in Joules) in the 
reference frame moving with the CM of the cloud. It 
follows from ([2"T)) that the rate of heating is proportional 
to Dp whose behavior is different in the regimes of regular 
and chaotic underlying Hamiltonian dynamics. 

The linear extent of the cloud in meters is Lx = 
2a x /kf. On a comparatively short time scale, r <C 
\((p c )/F)\, and low temperatures a v <C \(j?c)\, D p is ap- 
proximately the same for all the atoms in the cloud be- 
cause the CM velocity could not change significantly un- 
der the action of the friction force F during the observa- 
tion time. Using the first equation in the set ([3]) and the 
Eg. ([2"6"]l. we obtain 



1 



■ 2 4(0)r 2 



2,3 



-DpUJ^T 



(28) 



We have computed Lx with that formula with D p given 
by (|23[) and compare the result with numerical simula- 
tion of Eqs. ([3]) for a number of atomic clouds with dif- 
ferent initial values of (p c ). In Fig. [8] the dependence 
Lx((p c )) is plotted for A = —0.01 at two moments of 
time. The analytic dashed lines fit well the exact numer- 
ical results in the range (jp c ) < 1200 where the underlying 
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FIG. 8. A log-log dependence of the cloud's linear extent Lx 
(in microns) on the average initial momentum (p c ) of atoms 
in a cloud at two moments of time. The analytic dashed lines 
were computed with the formula 1)281) with D p = Dp valid 
in the chaotic regime on a short time scale. Note an abrupt 
change in the decay of Lx({Pc)) in the range of (p c ) ~ 1200 
where a chaos-order transition takes place in the underlying 
Hamiltonian motion. A = -0.01, 2a x (0) = 0.5, 2a p (0) = 5, 
wavelength A/ = 2n/kf = 850 nm. 



Hamiltonian dynamics is chaotic (see the left column in 
Fig. [7]). Note an abrupt change in the decay of Lx((p c )) 
beginning with those values of (p c ) ~ 1200 where the 
Hamiltonian dynamics becomes more regular. Since the 
linear extent of the atomic clouds changes abruptly at the 
chaos-order border one may conclude that in real experi- 
ments the value of Lx should increase sharply with those 
values of the average cloud momentum (p) for which the 
underlying Hamiltonian evolution is chaotic. 



VII. CONCLUSION 

Coherent evolution of the atomic state in a near- 
resonant standing-wave laser field is interrupted by SE 
events at random moments of times. The Hamiltonian 
evolution between these events has been shown previ- 
ously (for a summary of Hamiltonian theory for cold 
atoms in a ID optical lattice see Ref. [HI) to be reg- 
ular or chaotic depending on the values of the detuning 
A and the initial momentum po . We stress that dynam- 
ical chaos may happen without any noise and any mod- 
ulation of the lattice parameters. It is a specific kind 
of dynamical instability in the fundamental interaction 
between the matter and radiation. 

In reality Hamiltonian chaos is masked by random 
events of SE. The behavior of spontaneously emitting 
atoms in the detuning and momentum regimes where 
the underlying Hamiltonian dynamics is chaotic may be 
called stochastic chaos. We have specified and quanti- 
fied two regimes of the stochastic chaos, namely, random 
walking and dissipative ballistic transport. In the first 
regime, atoms wander in an optical lattice in a random 
way performing flights in both the directions with the 
PDFs strongly depending on the detuning (see Figs. [5] 



and [6]). In the ballistic regime, atoms move in the same 
direction with momentum fluctuations caused both by 
the Hamiltonian instability as well as SE events. It has 
been shown in our numerical experiments and confirmed 
analytically that the character of momentum diffusion 
changed abruptly in the regime where the underlying 
Hamiltonian dynamics proved to be chaotic. A clear cor- 
relation between the decay of the momentum diffusion 
coefficient D p and probability of Hamiltonian chaos A 
has been found (Fig. [7]). In order to observe the man- 
ifestation of Hamiltonian chaos in real experiments we 
proposed to measure a linear extent of atomic clouds Lx 
in a ID optical lattice and predicted a significant increase 
in Lx for the atomic clouds with A ~ 1. 

In conclusion wc would like to discuss some possible 
applications of the theory developed and the results ob- 
tained. A sensitive dependence of statistical properties of 
dissipative chaotic walking and ballistic transport on the 
values of the detuning A provides a possibility to manip- 
ulate atomic CM motion by changing A. For example, 
one can increase the mean duration of atomic flights in 
three orders of magnitude by changing A only by thirty 
percents (see Fig. |BJ). 

Cold atoms in optical lattices is an ideal system to 
study different phenomena in statistical physics. Besides 
dynamical chaos, the phenomena of stochastic resonance 
has been observed in a near-resonant optical lattice [46[ . 
Another phenomenon of considerable current interest is 
cold atom ratchets (47l - l5l| . A ratchet is a spatially peri- 
odic device which is able to produce a directed transport 
of particles in the absence of a net bias (i. e., when 
the time- and space-averaged forces are zero). In order 
to realize the ratchet effect it is necessary to break time 
or/and spatial symmetries which generate a countermov- 
ing partner to each trajectory [52(. Different classes of 
the ratchets have been experimentally realized with cold 
atoms in optical lattices [47l - l49j . The interrelation of 
Hamiltonian chaos and SE noise, found in this paper, 
provides an additional possibility to create and manipu- 
late directed transport of atoms in rigid optical lattices. 



APPENDIX 

We work in the regime of small detunings 

|A| < 1, (29) 
moderate mean atomic velocities 

(|w r p|) ~ 7 /2 « 1, (30) 
and diffusive motion 

r > 2/ 7 . (31) 



Due to (|29|) . we may neglect the last term in the poten- 
tial energy and adopt the Hamiltonian solutions for 
(UJ between any two acts of SE. The evolution of u is de- 
scribed by the approximate solutions ((5J) for the regular 
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Raman-Nath motion and J9|) for the chaotic motion. It 
follows from the condition of moderate atomic velocity 
([30]) and solutions (HI), © that |u| < 1. In other words, 
u cannot go far from zero between acts of SE after each 
of which u = 0. Under the conditions (|29|) and |u| <C 1, 
the mean kinetic atomic energy is much greater than the 
potential one and in this weak Raman-Nath approxima- 
tion 



In the chaotic regime we have from ([9]) the stochastic 
map 



ch = | A | 



■ sin ( 



UJ r p 



(36) 



Az cu r p 
\usmx + — <C — — 



the solution (|7|) for z is simplified 



z(t) ~ =p sin 



cos xdr' 



Xo 



(32) 



(33) 



In the diffusion regime (f3"Tj) and in the weak Raman-Nath 
approximation (|32p the momentum fluctuations between 
two SE are small in comparison with its fluctuations over 
the time scale of atomic transport. So, the atomic mo- 
mentum just before SE is equal to its value at the node 



(34) 



Now we can get simplified solutions for u. In the ex- 
act Raman-Nath condition, x = uj r pr, we have from ([5]) 
the approximate deterministic map written in a two-step 
form 



,RN 



2A 



uj r p 



Vq cos 



2 

UJ r p 



,RN 



(35) 



The maps derived enable us to estimate the values of 
u(tJ) just before jth act of SE. After the (j - l)th SE 

uq = u(rt x ) = 0, u(t~) = um is an accumulated value 
of u after passing M nodes in the interval Tj_i < r < tj. 
The average number of node crossings can be estimated 
to be (M) ~ 2\oj r p\/jTT. In the exact Raman-Nath limit, 
um = M(u m — u m _2)/2, and the mean squared value in 
the regular regime 



{ u2 { T J)) RN = «i) 



L MIRN 



4A 2 clvp , 9 . o/2 

5— (Vq) cos z 

7T7 2 u/ \u} r p 4 

(37) 

Beyond the strong Raman-Nath limit, u[t^ ) = um is 
a sum of M random numbers which are proportional to 
sin (j} m . From the probability theory we get 



2u! r p 



A 2 

— (38) 
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